CMB polarization features from inflation versus reionization 
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The angular power spectrum of the cosmic microwave background temperature anisotropy ob- 
served by WMAP has an anomalous dip at ^ ~ 20 and bump at ^ ~ 40. One explanation for this 
structure is the presence of features in the primordial curvature power spectrum, possibly caused 
by a step in the inflationary potential. The detection of these features is only marginally significant 
from temperature data alone. However, the inflationary feature hypothesis predicts a specific shape 
for the iS-mode polarization power spectrum with a structure similar to that observed in temper- 
ature at ^ ~ 20 — 40. Measurement of the CMB polarization on few-degree scales can therefore 
be used as a consistency check of the hypothesis. The Planck satellite has the statistical sensitiv- 
ity to confirm or rule out the model that best fits the temperature features with 3 a significance, 
assuming all other parameters are known. With a cosmic variance limited experiment, this signifi- 
cance improves to 8 a. For tests of inflationary models that can explain both the dip and bump in 
temperature, the primary source of uncertainty is confusion with polarization features created by a 
complex reionization history, which at most reduces the significance to 2.5 a for Planck and 5 — 6 cr 
for an ideal experiment. Smoothing of the polarization spectrum by a large tensor component only 
slightly reduces the ability of polarization to test for inflationary features, as does requiring that po- 
larization is consistent with the observed temperature spectrum given the expected low level of TE 
correlation on few-degree scales. If polarized foregrounds can be adequately subtracted, Planck will 
supply valuable evidence for or against features in the primordial power spectrum. A future high- 
sensitivity polarization satellite would enable a decisive test of the feature hypothesis and provide 
complementary information about the shape of a possible step in the inflationary potential. 



I. INTRODUCTION 



Our best constraints on the shape of the primordial 
power spectrum at large scales come from observations 
of the cosmic microwave background (CMB) anisotropy 
by the Wilkinson Microwave Anisotropy Probe (WMAP) 
HI, [2. The latest (5-year) WMAP data continue to 
be well described by the simplest inflationary scenario 
of a single, slowly rolling, minimally coupled scalar field 
with a canonical kinetic term @, 0, H, Q ■ Since the 3- 
year release [l3| , the WMAP data have indicated a devi- 
ation from scale invariance — a red tilt of the scalar spec- 
tral index — the significance of which has been debated 
in the literature from a Bayesian model selection point of 
view (e.g. (Til, fl^). Recent minimally-parametric recon- 
structions of the primordial power spectrum incorporat- 
ing some form of penalty for "unnecessary" complexity 
[l3l . [T3 | show some evidence for a red tilt, but no evi- 
dence for scale dependence of the spectral index. These 
methods, as currently implemented, are not very sensi- 
tive to sharp, localized features in the primordial power 
spectrum. 
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However, it has been pointed out ever since the orig- 
inal data release [1^, [lB| that there arc several sharp 
glitches in the WMAP temperature (TT) power spec- 
trum. In particular, several model-independent recon- 
struction techniques that are sensitive to features local- 
ized in a narrow wavenumber range have consistently 
picked out a feature at £ ~ 20 — 40 that leads to an 
improvement of ~ 0(10) over a smooth power-law 
spectrum [13, M, \mM,^- 

Power spectrum features could arise, in principle, in 
more general classes of inflationary models where slow 
roll is momentarily violated. Such an effect can be phe- 
nomenologically modeled as a discontinuity or singularity 
in the inflaton potential [H, [H, [111 . A "step- like" fea- 
ture [1^, in particular, would be a good effective field 
theory description of a symmetry breaking phase transi- 
tion in a field coupled to the inflaton in multi-field models 
m, m, 113, m, HI, which can arise in supergravity [s^l 
or M-theory- inspired [Sll, [13 contexts. Several analyses 
have confronted such phenomenological descriptions of 
features in the inflationary potential with current data 

[H, M M MMM,Mm- 

It is debatable whether the large scale feature seen in 
the WMAP TT spectrum is a signal of exotic primor- 
dial physics or merely a statistical anomaly. Currently, 
our information about the smoothness of the primordial 
power spectrum is dominated by the temperature data. 
However, future high fidelity CMB polarization measure- 
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mcnts at large scales have the potential to shed light on 
this question. The importance of polarization data for 
constraining oscillatory features has be en p reviously dis- 
cussed in the literature (e.g. [Ill, [13, SE S IH) and 
exploited in particular as a cross-check of the observed 
low CMB temperature quadrupole [13, lil, lil, [i^l • 

In this work, we propose to use the large-scale po- 
larization of the CMB to test the hypothesis that the 
I ^ 20 — 40 glitch is due to a step in the inflaton po- 
tential. We exploit the fact that, in the relevant mul- 
tipole range, the sharpness of the polarization transfer 
function and lack of contamination by secondary effects 
(assuming instantaneous reionization) makes polariza- 
tion a cleaner probe of such features than temperature 
j40| . We also investigate how our conclusions are affected 
by rela xing the assumption of instantaneous reionization 
psi . [49I |50| | . changing the parameters of the feature, and 
including large- amplitude tensor fluctuations. This anal- 
ysis is particularly timely given the imminent launch of 
the Planck satellite [5l|, which promises to greatly in- 
crease our knowledge of the large-scale polarization sig- 
nal. It is also relevant for future dedicated CMB polar- 
ization missions [52j . As in a related previous paper on 
polarization consistency tests of large-angle CMB tem- 
perature anomalies [53| . it is our objective to make a 
prediction for the polarization statistics that will be ob- 
served by future CMB experiments, given current tem- 
perature data, in a "last stand" before Planck. 

We present the inflationary model and the numeri- 
cal procedure used to compute the primordial curvature 
power spectrum in § [III The polarization consistency 
tests of the features, both for instantaneous and general 
reionization histories, are presented in § IIIII and § IIVI 
and we conclude in § [V] We discuss in Appendix [X] the 
relation of our work to previous analyses of features in 
the WMAP temperature data. 



with the amplitude and width of the step determined by c 
and d respectively, assuming that both are positive num- 
bers. We express the potential parameters m, 6, and d in 
units of the reduced Planck mass, Afpi = (SttG)^^/^ = 
2.435 X 10^^ GeV; the step amplitude c is dimcnsionlcss. 

In physically realistic models with a sufficiently small 
step in the potential, the interruption of slow roll as the 
field encounters the step docs not end inflation but affects 
density perturbations through the generation of scale- 
dependent oscillations that eventually die away. The 
phenomenology of these oscillations is described in Ref. 
[23|: the sharper the step, the larger the amplitude and 
width of the "ringing" superimposed upon the underly- 
ing smooth power spectrum. Hence we shall see in ij llllBI 
that lowering d increases the width of the feature in I in 
the CMB power spectra. 

Standard slow-roll based approaches are insufficient for 
computing the power spectrum for this potential, and 
instead the equation of motion must be integrated nu- 
merically mode- by- mode js^l- It is convenient to use the 
gauge invariant Mukhanov potential [55l.[56| for the mode 
amplitude since it is simply related to the curvature per- 
turbation IZ: 

u ^ -z7^ , (2) 

where z = 4>/H, H is the Hubble parameter, and the 
dot denotes a derivative with respect to conformal time. 
The Fourier components Uk obey the equation of motion 

[13, 

Uk + (^fc' - ^ j ufc - , (3) 

where k is the modulus of the wavevector k. The power 
spectrum is defined via the two point correlation function 



II. INFLATIONARY FEATURES 



We review the inflationary generation of features in the 
curvature power spectrum from step-like features in the 
inflaton potential in miA\ and their transfer to the CMB 
temperature and power spectra in ijllBI 



which is related to Uk and z via 



(5) 



A. Inflationary Model 

To model a feature in the primordial power spectrum 
that matches the glitches in the WMAP temperature 
data at ^ ^ 20 — 40, we adopt a phenomenological infla- 



tionary potential of the form V^(0) = m^fj((/))(/) /2 where 
the effective mass of the inflaton (p has a step at cf) = b 
corresponding to the sudden change in mass during a 
phase transition [2^: 



'eff 



(0) 



;tanh 



(1) 



The dynamics of the Hubble parameter, described by 
the Friedmann equation, and the background dynamics 
of the unperturbed inflaton field, described by the Klein- 
Gordon equation, can be written respectively as 



H' 



H' 



1 dV 



= 0, 



(6) 
(7) 



where ' = d/d\na. The solution of the mode equation 
depends on the background dynamics. With the help of 
these background equations, the mode equation ([3]) can 
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TABLE I: Fiducial feature model parameters chosen to best 
fit WMAP5 under a flat ACDM cosmology, compared in the 
text with a smooth model with c = 0, m = 7.120 x 10~*, and 
the same cosmological parameters, which matches the small 
scale normalization As(fc*) — 2.137 x 10~® and tilt Us ~ 0.96 
at the pivot fc* = 0.05 Mpc"\ 



Parameter 


Value 


m 


7.126 X 10"^ 


6 


14.668 


c 


1.505 X 10"^ 


d 


0.02705 




50 




0.02238 




0.1081 


h 


0.724 


T 


0.089 



FIG. 1: Upper panel: solid black: Inflationary potential with a step 
[Eq. [T]. The parameters for the potential are chosen to maximize 
the WMAP5 likelihood and arc listed in Table |l] The dashed red 
line shows a smooth m?(j>^ potential (c = 0) with m = 7.120 X 10^'' 
so that the two models have equal power on small scales (4> ^b). 
Middle and lower panels: slow-roll parameters ey and r)v for the 
two inflationary potentials. 



be written as 




(8) 



where the term in square brackets is z/(za?'H'^). We stop 
integrating the background equations after any transient 
solution has died away, but while the mode is still well 
within the horizon. This allows us to obtain initial con- 
ditions for the two orthogonal solutions that contribute 
to Uk that are free of contamination due to any transient 
contribution to the background dynamics. The power 
spectrum is then obtained by continuing the integration 
until the mode freezes out far outside the horizon, yield- 
ing the asymptotic value of Further details re- 
garding the numerical solution of the coupled system of 
differential equations can be found in Ref. [23j . 

To match a given mode to a physical wavenumber 
k, one must make an assumption about the reheating 
temperature, but this choice is degenerate with &, corre- 
sponding to a translation of the s tep in 0. To compare 
our results with those of Refs. [H, [3J|, we adopt the fol- 
lowing prescription for the matching: 



where H^, is the Hubble scale corresponding to the phys- 
ical wavenumber A:^,, which left the horizon N^, e-folds 
before the end of inflation, defined by (Pa/dt^{acnd) = 0. 
Following the above authors, we set the pivot scale 
k^, = 0.05 Mpc~^ to correspond to A^^ = 50 (although 
there are differences in the implementation of the fc-mode 
matching that we discuss in Appendix VK^ . 

Figure [T] shows our fiducial inflationary potential, with 
parameters given in Table H] that are chosen to fit the 
WMAP5 temperature glitches at ^ ~ 20 — 40 as we will 
show in the next section. The number of e-folds of infla- 
tion after the step in this potential is A^step ~ 54. The 
slow-roll parameters 



PI (dV/d(t>\ 



Ml, 



(PV/df_ 
V 



(10) 



(9) 



are plotted in the lower panels of Fig. [T] Note that near 
the step at = 6, |77y| > 1 confirming that the slow-roll 
approximation is not valid. Figure [5] shows the inflation- 
ary curvature power spectrum A^(fc) for this potential, 
computed by integrating Eqs. ([SI)—®. 

For comparison, in Figs. [1] and [5] we also show a 
smooth, c = potential with the same small-scale ampli- 
tude and tilt as the fiducial potential, and its slow-roll pa- 
rameters and inflationary power spectrum. The smooth 
spectrum is nearly indistinguishable from a pure power 
law of Us ~ 0.96 with amplitude As(fc^) = 2.137 x 10"^. 
Note that the spectral index is determined by the choices 
of Ni, and fc^ in the matching condition of Eq. ([9]), while 
the amplitude comes from the inflaton mass m. 
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FIG. 2: Primordial curvature power spectra for the potentials in 
Fig.E 

B. CMB Power Spectra 

The mapping between the inflationary curvature power 
spectrum and the observable CMB angular power spectra 

{^em^e'm') = ^U'^mm'Cf^ , (11) 

where X, X' e T, E, is given by the scalar radiation 
transfer functions 

^i^±il^ = |rflnfcTf(fc)T/'(fc)A^(fc). (12) 

In Fig. [21 we show the T and E transfer functions for 
the fiducial cosmological parameters of Table HI For a 
more extended discussion of the transfer functions and 
their relationship to features in the inflationary power 
spectrum, see j4C| . The resultant temperature and polar- 
ization angular power spectra from the inflationary power 
spectra of Fig. [2] are plotted in Fig. [H 

For the wavenumbers of interest, 1 < fc/lO^'^Mpc"^ < 
5, the transfer of power to temperature fluctuations tran- 
sitions between the Sachs- Wolfe and acoustic regimes at 
high £ and carries substantial contributions from the in- 
tegrated Sachs- Wolfe (ISW) effect at low £. These effects 
and geometric projection lead to a very broad mapping 
of power in k to power in i. In particular, the oscillations 
at the upper range in k are largely washed out, leaving 
only a single broad dip at ^ ~ 20 and bump at £ ~ 40 in 
the temperature spectrum. Likewise, the power at these 
multipoles correspond to a wide range in k as shown in 
Fig.El 

Polarization spectra differ notably from the tempera- 
ture spectra due to the differences in the transfer function 
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FIG. 3: Transfer function (k) for the fiducial model with in- 
stantaneous reionization. Upper panel: temperature X = T; lower 
panel: polarization X = E. Contours are spaced by factors of 2. 
Dashed lines represent the range of fe-modes where features appear 
in Fig. [2] Polarization is a cleaner probe of features in this range 
and, for instantaneous reionization, is nearly uncontaminated by 
secondary effects. The temperature and polarization are also only 
weakly correlated here due to the transition between the Sachs- 
Wolfe (SW) and acoustic regimes in temperature. 



shown in Fig. [31 For the standard instantaneous reion- 
ization history and the upper portion of the range of k 
affected by the feature, the polarization is dominated by 
the onset of acoustic effects only. We shall see that this 
makes the bump in ~ 40 a particularly clean test of in- 
flationary features (see Fig. [S]) . Furthermore oscillations 
from high k at higher £ are retained at a significant level 
in the polarization. 

On the other hand at fc ^ 10"'^ Mpc~^, the polariza- 
tion transfer from recombination becomes very inefficient 
and reionization effects come into play. This leads to a 
very low level of polarization around £ ~ 20 with features 
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FIG. 4: Temperature and polarization power spectra for the in- 
flationary power spectra in Fig. [2] with solid black curves for the 
model with a feature and dashed red curves for smooth A'^{k). 
Dotted curves indicate where Cj^ is negative. Blue points with 
error bars show the 5-year WMAP measurements of Cf^ includ- 
ing sample variance. For both models, the reionization history is 
assumed to be instantaneous and the cosmological parameters not 
determined by the inflationary potential are given in Table |l] 



even for a smooth inflationary power spectrum. These 
properties leave the ^ ~ 20 dip vuhierable to external 
contamination such as tensor contributions (see § IIII C[) 
or foregrounds as well as uncertainties in the ionization 
history (see § lIVp . 

Finally, the cross correlation between the temperature 
and polarization fields for the entire range of 20 < £ < 40 
is very low due to the transition between the Sachs- Wolfe 
and acoustic-dominated regimes in the temperature field. 
We shall see in ij llll DI that this prevents statistical fluctu- 
ations in the observed temperature power spectrum from 
being repeated in the polarization. 

These differences in the transfer functions also play 
a role in defining the region in the potential parameter 
space that best fits the WMAP TT data versus the re- 
gion that is best tested by polarization. For the former, 
we conduct a grid based search over the potential param- 
eters. The mass parameter m determines the amplitude 
of the spectrum away from the feature and so is mainly 
fixed by the acoustic peaks at high £. The location of the 
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FIG. 5: Transfer function (k) for the fiducial model with in- 
stantaneous reionization for multipoles near the temperature dip 
{1 = 20) and bump (£ = 40) for temperature and polarization. For 
temperature, the dip multipoles receive a broad range of contribu- 
tions from k > 10~^ and the bump multipoles from A; > 3 X 10^'^. 
The localization of the transfer function is sharper for polarization, 
especially for £ = AO which is immune to reionization effects. The 
polarization transfer functions have been scaled by IC* and 10''' for 
convenience. 

feature b is also well determined independently of the 
other parameters [H, [U . We therefore fix m and b at 
their best-fit values and search for the best fit in the step 
amplitude and width parameters c and d. 

The values of to, b, c, and d given in Table H] specify 
the maximum likelihood model. This model improves 
the fit to the 5-year WMAP data by -2A In ^tt « -8. 
We will explore variations in the parameters about the 
maximum and their relationship to the temperature and 
polarization power spectra through the transfer functions 

The improvement is only marginally significant given 
the 3 extra parameters of the step and the choice of one 
out of many possible forms. Matching polarization fea- 
tures can therefore provide a critical confirmation or refu- 
tation of the inflationary nature of the temperature fea- 
tures. 



III. CONFIRMING FEATURES WITH 
POLARIZATION 

In this section, we discuss the significance with which 
polarization measurements can confirm or rule out the in- 
flationary features discussed in the previous section under 
the instantaneous reionization model. We begin in S IIII Al 
with the significance of the best-fit feature model under 
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TABLE II: Parameters used when making forecasts for ide- 
alized and Planck-like experiments. Here Ap' is in units of 
^K-arcmin. 



Experiment 


V 


■^FWHM 




/sky 


Ideal 










0.8 


Planck 


70 GHz 


14.0' 


255.6 


0.8 




100 GHz 


10.0' 


109.0 


0.8 




143 GHz 


7.1' 


81.3 


0.8 



TABLE III: — 2Aln£_B_E for false positive and false negative 
tests comparing models with smooth A|j^(fc) and a feature in 



Experiment Test 




Ideal 


False positive 


64 


Ideal 


False negative 


60 


Planck 


False positive 


8 


Planck 


False negative 


9 



the simplest set of assumptions. We assess changes in the 
significance due to variation in the potential parameters 
in § IIII Bl and due to the inclusion of tensor i?-modcs 
in § IIII CI In § IIIIDl we describe the impact of condi- 
tioning polarization predictions on the already-measured 
temperature spectrum. 



Table Hill lists —2^\x\Lee for rejecting false posi- 
tives and false negatives. The significance of false posi- 
tive or negative rejection in this most optimistic case is 
2A \wLee ^ 8 for the ideal experiment and ~ 3 for 
Planck. In the following sections, we will discuss various 
effects that can degrade this significance. 



A. Fiducial Polarization Significance 



B. Potential Parameters 



To evaluate the significance of discriminating between 
models, we assume a Gaussian likelihood for the polariza- 
tion angular power spectrum. In the absence of detector 
noise, the likelihood Lee of data Cf^ given a model 
power spectrum Cf^ is 



IX'uLee 



1) 



EE 



EE 



In 



C, 



EE 



EE 



- 1 



(13) 

where /sky is the fraction of sky with usable E measure- 
ments. If the data Cf^ have no inflationary feature and 
the model Cf^ has the inflationary feature, we call this 
the significance at which false positives can be rejected. 
Conversely, if the data have an inflationary feature and 
the model spectrum has no feature, we call this the sig- 
nificance at which false negatives can be rejected. 

For forecasts throughout this paper, we assume that 
the data are equal to the ensemble average of realiza- 
tions for a particular model. Therefore, the minimum 
— 21n£ is zero and — 2Aln£ = — 21n£. The exception 
to this is that when we discuss the WMAP TT likelihood, 
the relevant quantity is — 2Aln£ = — 21n(£/£ML) where 
the likelihood of the best fit model is — 21n£ML 7^ 0. 
We make forecasts for an ideal, sample variance limited 
experiment and for Planck using the experimental spec- 
ifications in Table [III For the Planck case with a finite 
noise power N^, Ce ^ Ce + in Eq. where Ng is 

the minimum variance combination of the noise powers 
of the individual frequency channels 



A 



^K-rad 



exp 



l)(4wM/^ad)^ 
81n2 



(14) 



using Ap and ^'fwhm from Table |TT] converted to the 
appropriate units. 



Variation in the parameters of the inflationary poten- 
tial from the best flt model can affect the significance of 
polarization tests of features. As noted in § III B[ m and 
b are strongly constrained by the observed CMB temper- 
ature spectrum, but the parameters c and d that control 
the amplitude and width of an infiationary step are less 
well determined by temperature alone. 

In terms of the curvature power spectrum, increasing 
c increases the amplitude of the features. However, de- 
creasing the width of the potential step by lowering d 
enhances the deviations from slow roll, thereby also am- 
plifying the feature in the power spectrum. 

Figure [5] shows a contour plot of the WMAP temper- 
ature likelihood —2AlnCTT for the parameters c and d 




0.001 



0.002 



0.003 



FIG. 6: Contour plot of — 2A In Ctt for parameters c and d using 
5-year WMAP data. Other potential parameters are fixed at their 
fiducial values. The minimum, with —2Aln Ctt = —8.3 relative 
to the smooth c = model, is shown with a cross. 
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FIG. 7; Contour plot of —2Aln Cee for the parameters c and d, 
for tests of false positives with a cosmic variance limited experiment 
(upper panel) and Planck {lower panel). The best fit model to 
WMAP TT is shown with a cross. 

(relative to c = 0) and Fig. [7] shows — 2 A\n Cee for false 
positives using simulated polarization data. The simi- 
larities and differences between these two plots reflect 
properties of the temperature and polarization transfer 
functions. 

For the temperature case near the minimum, the de- 
generacy between the two parameters is approximately 
c ex. (P. This line roughly corresponds to keeping the 
amplitude of the enhanced power in A^(fc) at fc ^ 
3 X 10"'^ Mpc~^ fixed. The preferred value of c cor- 
responds to the best amplitude of the negative dip at 
k ^ 2 X 10~^ Mpc~^. For the best fit parameters, in- 
cluding the feature in A|j(fc) improves the fit to 5-ycar 
WMAP data by -2 A In ^tt ~ -8. 

Due to the weak significance of the feature detection, 
the contours become substantially distorted away from 
the maximum likelihood. In particular, the contours in 
Fig. [6] show a triangular region extending to high d ^ 
0.04. This region corresponds to a lower amplitude in 
both the first dip and bump in k as shown in Fig. [H] Due 
to projection effects in temperature, the £ = 20 dip gets 
contributions from both the dip and the bump in k (see 
Fig. [51) . Consequently, a model with smaller features in k 
in both the dip and bump can lead to the same amplitude 
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FIG. 8: Primordial curvature power spectra for models illustrating 
projection degeneracies in the temperature. The parameters of the 
two models are chosen to have similar temperature dips at £ ~ 
20 and equal WMAP likelihoods: {c,d) = (0.00128,0.043) (solid 
black) and (c,d) = (0.0023,0.028) (dashed red); other parameters 
are fixed to the values in Table |T] for both models. 
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FIG. 9: Upper panel: Relative difference in Cj"^ , with respect to 
a smooth power spectrum, of two models with equal TT amplitude 
at £ ~ 20 and curvature power spectra shown in Fig. [8] Lower 
panel: Due to differences in the polarization transfer function, the 
models do not have degenerate £ ~ 20 dips and show significant 
differences at £ > 60 as well, leading to a much higher polajization 
significance for the model plotted with dashed lines. 



of the dip at £ = 20 if the amplitude of the bump is 
reduced more. 

We illustrate these projection effects in Figs. [5] and [5] 
with two models chosen to have the same likelihood im- 
provement of — 2Aln£TT w —6. For the model with 
smaller features in k, the temperature enhancement at 
^ ^ 40 is substantially reduced compared with the best- 
fit model, while the model with larger features in k over- 
shoots the bump at £ = 40 in temperature. Despite these 
differences, both models have about the same TT ampli- 
tude in the £ = 20 dip as the best-fit model but a slightly 
worse overall fit. In particular, for the model with smaller 
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features in k, inflationary features can only explain the 
observed £ = 20 dip in temperature and not the £ = 40 
bump. 

The degeneracies in c and d for polarization signif- 
icance share similarities with, yet have important dif- 
ferences from, those for temperature. The polarization 
significance remains largely unchanged for small varia- 
tions in c and d along the constant c/d^ line favored 
by the temperature spectra. Near the maximum, vari- 
ations along this direction preserve the amplitude of in- 
trinsic features in k (see Fig. [5]). However within the 
— 2Aln£7-7- = —4 region the significance for a ideal ex- 
periment can either drop or rise significantly. The rea- 
son is that due to projection effects in temperature, the 
polarization better separates changes in the overall and 
relative amplitude of the features in fc. In the triangular 
high d region, where the amplitude of the TT dip remains 
unchanged but the intrinsic features in k are all reduced, 
the significance of the polarization difference decreases 
markedly (see Fig.lHl). Because of the sharper projection, 
even the £ ~ 20 dip in polarization is reduced. The net 
result is that the polarization significance is a stronger 
function of c, which controls the overall amplitude, than 
the temperature significance. 

Note that while the significance can be substantially 
degraded from our best fit assumptions, this is mainly 
because of the weak detection of a feature in the tem- 
perature spectrum itself. In cases where the polarization 
significance is greatly reduced, the temperature bump at 
£ ~ 40 cannot be explained by the infiationary features. 
In other words, polarization remains a robust probe of 
the infiationary nature of the £ = 40 bump across varia- 
tions of the potential parameters. 



C. Tensors 

The m^0^ potential with the parameters in Tablc|T]pre- 
dicts substantial gravitational wave contributions with 
tensor-to-scalar ratio r w 0.16. Relative to a smooth- 
A^(A:) model without tensors, the to^0^ model with a 
feature has extra distinguishing power due to the pres- 
ence of i3-mode polarization. Because other forms for the 
potential can also be used as the smooth base on which 
to place the feature [11] we choose not to include tensors 
for most of our calculations. Moreover, a B-mode de- 
tection would not be useful for discriminating features. 
The _B-mode amplitude is insensitive to features since the 
potential amplitude is left nearly unchanged by the step. 
Additionally, a small step in the potential is not expected 
to generate features in the CMB power spectra of tensor 
modes. Unlike the scalar spectrum whose shape is sen- 
sitive to the second slow- roll parameter r]v, the shape 
of the tensor spectrum depends primarily on ey, which 
remains small at the step (see Fig. [U (S^). In Fig. [TU] 
we show the i?-mode prediction for r = 0.16 and a pure 
power law tensor spectrum with tilt nt = — r/8. 

On the other hand, it is important to assess the possi- 
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FIG. 10: Effect of tensor fluctuations on polarization power spec- 
tra for the model with a feature in A^(fc). Solid black: no tensor 
component. Dashed red: including tensors with r = 0.16. Tensors 
smooth the EE spectrum near the £ ^ 20 dip. 



bility of degradation of the E-mode feature from the cur- 
vature spectrum due to the nearly smooth tensor i?-mode 
contributions. Due to the shape of the tensor i?-mode 
spectrum, which mimics the B-mode spectrum, the main 
impact of tensors is to fill in the dip in the polarization 
spectrum around £ ^ 20 (sec Fig. [T0|. Correspondingly, 
the decrease in significance for the ideal experiment is 
13 - 15% in y/-2A\nCEE, and for Planck, 4%. Planck 
is less affected since its lower sensitivity limits the ac- 
curacy of measurements in the £ ~ 20 dip. Since these 
degradations are relatively small, we ignore tensors when 
considering the impact of the reionization history below. 
Moreover, we shall see that reionization uncertainties are 
very similar to tensors in that they make the £ 20 dip 
less useful for distinguishing features through _B-mode 
polarization. 



D. Temperature Conditioning 

The usefulness of polarization for providing an inde- 
pendent test of features observed in temperature may 
also be reduced by the correlation of temperature and po- 
larization: a positive correlation would make observation 
of polarization features more likely given the WMAP TT 
data regardless of whether the features have an inflation- 
ary or chance statistical origin. We expect the reduction 
in significance to be small given that Cf^ is small on the 
relevant angular scales (see Fig. U]) , and in this section 
we quantify this statement. 
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To assess the impact of conditioning polarization pre- 
dictions on the WMAP temperature data, it is convenient 
to replace the likelihood statistic of Eq. (fT3|) with a ^ 
statistic. This allows us to phrase the impact in terms 
of the bias and change in variance predicted for the EE 
power spectrum from the TT measurements. Note that 
in the absence of the temperature constraint and in the 
limit of small differences between the model and the data, 
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which is equal to a simple statistic. 

Now let us include the temperature constraint. First 
take the idealization that the temperature multipolc mo- 
ments Tern have been measured on the full sky with neg- 
ligible noise. Given a model that correlates the polariza- 
tion field through the cross correlation coefficient 
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a constrained realization of the polarization field that is 
consistent with the temperature field can be constructed 
as 



Re 



EE 



c\ 



TT 



1 - -Rf fffn 



(17) 



where gg^a is a complex Gaussian field with zero mean, 
unit variance {gem9grn) = 1' ^^"^ ^ ^'^^^ transform g^^^ = 
(~l)™(7^._„i. The estimate of the power spectrum is then 
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and its mean over the constrained realizations is biased 
from the true 
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by the fixed observed temperature power spectrum 
CF = Y.mTtmTUi'^^ + !)■ With a high correlation 
coefficient, chance features in the temperature spectrum 
induce similar features in the observed polarization spec- 
trum. For example, if TT fluctuates high, EE will also 
fluctuate high (on average). 

The temperature constraint also removes some of the 
freedom in the variance of the polarization power spec- 
trum: 
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FIG. 11: Upper panel: Solid lines show the _E-mode power spec- 
trum constrained to the temperature data for the smooth- A|j(fc) 
model along with the band representing sample variance per I for 
the ideal experiment. The model with a feature {dashed) lies signif- 
icantly outside of the band in the 10 < ^ < 60 range, making false 
negatives unlikely. Lower panel: Fractional difference between the 
average of the constrained realizations {Cf^) and the full ensem- 
ble average for both models. The impact of the constraint 
is minimal due to the lack of correlation between the temperature 
and polarization fields in the region of interest. 



In the limit that the correlation R^ — 0, the variance 
takes on its usual form for a Gaussian random field. In 
the limit that i?f ^ 1, there is no uncorrelated piece 
and the observed temperature spectrum determines the 
observed polarization spectrum with no variance. 

Now let us add in detector noise and finite sky cover- 
age. Given a noise power spectrum f-^d a fraction 
of the sky /^ky. 
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FigurefTTlshows. in the upper panel, the i?-mode polar- 
ization power spectrum for the smooth inflationary spec- 
trum constrained to WMAP5 temperature data for the 
ideal experiment. For comparison, Cf^ for the best fit 
feature model is also plotted. Note that even the sec- 
ond dip in the spectrum at £ ^ 60 remains significantly 
distinct in polarization. In the lower panel, the impact 
of the temperature power spectrum constraint is plotted 
as the fractional difference between Cf^ and Cf^ for 
each model. Due to the lack of temperature-polarization 
correlation in the 10 < £ < 60 regime, the impact of the 
constraint on the polarization features is negligible. 

We can quantify these conclusions by generalizing the 
statistic in Eq. (jlSp to include the temperature con- 
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TABLE IV: Axbb for false positive and false negative tests 
comparing models with smooth A^lk) and a feature in 
A^(A;), with polarization either unconstrained or constrained 
to observed temperature data. 



^Xee 

Experiment Test w/o T with T 

Ideal False positive 70 63 

Ideal False negative 59 56 

Planck False positive 8 8 

Planck False negative 9 8 



TABLE V: — 2Aln£_B_B for tests of false positives and false 
negatives with ionization histories of the data and model 
tuned at 6 < z < 50 to minimize the significance of rejec- 
tion using the methods described in § IIVBI 



Experiment Test 


-2A\n Cee 


Ideal 


False positive 


36 


Ideal 


False negative 


25 


Planck 


False positive 


6 


Planck 


False negative 


6 



straint: 

AyIp = V — . (22) 

As in the likelihood analysis, we assume that the data are 
a typical draw of the true model ("1") and that we are 
testing the significance at which the second model ( "2" ) 

can be rejected. Then we set Cf ^ = (Cf ^^^^), (Cf ^) = 
(Cf^^'^), and Var(Cf^)= Var(C'f ^^''). Note that the 
bias induced by the temperature constraint enters into 
both models whereas the change in the variance enters 
only from model 2. 

Table HVl assesses the significance of the rejection of 
false positives and false negatives for the fiducial feature 
model. In the last column we have applied the constraint 
from the WMAP5 temperature data and in the penulti- 
mate column we artificially drop the constraint by setting 
Ri = in the evaluation of Axee- Even with the con- 
straint, the significance with which false positives can be 
rejected is ^yAxj^ = 7.9 for the ideal experiment and 
-\/Ax|7g = 2.8 for Planck. For the case of false negatives, 
these numbers become 7.5 for the ideal experiment and 
2.8 for Planck. The difference between false positive and 
false negative significances comes from the dependence of 
sample variance on the model tested. In all cases, the sig- 
nificance in terms of Axee comparable to — 2 Aln Cee 
in Table Uni 

The impact of the temperature constraint is to lower 
the significance of both cases but only by < 5% in 
y^^Ax^- This small difference in significance justifies 
our choice to omit the constraint to temperature in our 
exploration of other effects that can degrade the signifi- 
cance. 



IV. REIONIZATION FEATURES 

A more complicated ionization history can in principle 
produce features in the polarization spectrum that might 
mimic or obscure features from the inflationary power 
spectrum. This is especially true for the dip at ^ ^ 20. 



In this section we search for ionization histories that lead 
to a higher incidence of false positives and false negatives. 

Reduced significance of false positives or negatives due 
to confusion between features from inflation and from 
reionization can arise in two ways. First, the true reion- 
ization history can introduce features in the data that 
either falsely mimic inflationary features or hide true 
features. Second, additional reionization freedom in the 
(false) model we wish to test can allow a better match 
to data generated from the alternate (true) model as- 
sumption. To account for both effects, we use a two-step 
method in which we first optimize the ionization history 
of the true model to produce a false positive or negative 
result, and then vary the ionization history of the false 
model. We will describe this procedure in § IIVBI 

Table |V] summarizes the results of this section. Rela- 
tive to the significance of rejecting false positives or neg- 
atives for instantaneous reionization models, ionization 
freedom lowers the significance for an ideal experiment 
by a factor of 0.64 — 0.75 to ^/—2A IuCee ~ 5 — 6, and 
for Planck by a substantially smaller factor of 0.83 — 0.87 
to \/—2A\nCEE ~ 2.5. Given the amount of freedom 
we allow in the ionization history these should be viewed 
as the maximal degradation possible due to reionization. 
We describe the details of this calculation in the following 
sections. 



A. Reionization Principal Components 

The form of the ionization history, and therefore the 
shape of the large-scale reionization peak in the polar- 
ization spectrum, are only weakly constrained by current 
observations and theoretical modeling, especially on the 
scales relevant for inflationary features [601 ■ We treat 
the evolution of the mean ionized fraction of hydrogen 
with redshift, Xe{z), as an unconstrained function be- 
tween z = 6 and some high redshift z = Zmax- At lower 
redshifts, we assume « 1 as required by the observed 
Lyc t transmission in quasar spectra at z < 6 (see e.g. 
p3l|). The highest redshift of reionization is less certain, 
so we take Zmax = 50 which is quite conservative for con- 
ventional sources of ionizing radiation. 

We parametrize general reionization histories with a 
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basis of principal components (PCs) Si{z) of the large- 
scale E-mode polarization We use the 7 lowest- 
variance PCs only since the higher-variance PCs have 
a negligible impact on the polarization power spectrum. 
Thus the ionization history at 6 < z < 50 is 



i=l 



(23) 



of 




where we take a constant fiducial ionized fraction 
x^'^ = 0.07 so that the fiducial model with {m^} = 
has a total reionization optical depth of r w 0.09. 

We vary the PC amplitudes {m^} and compare the 
resulting CMB power spectra with data simulated for 
ideal and Planck experiments using Markov Chain Monte 
Carlo (MCMC) likelihood analysis as we describe in the 
next section. 



B. Data and Model Optimization 
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We use a two-step optimization process to determine 
the maximum reduction in significance of false positive 
or negative rejection that can be caused by reionization 
features. We categorize models in this section by whether 
A^(fc) is smooth (S) or has a feature (F), and by whether 
the reionization history is instantaneous (I) or more com- 
plex (C) and parametrized by principal components as in 
Eq. (p3)) . For comparisons of models we introduce the no- 
tation f alse:true; for example, FLSI represents the false 
positive test for instantaneous reionization models from 

In the case of false positives, the goal of the optimiza- 
tion is to go from the FLSI comparison of the previous 
section to FC:SC, in which both false and true models 
have complex ionization histories. In particular, we want 
to find the FC:SC pair that minimizes the difference be- 
tween the two models, thus minimizing the significance 
of false positive rejection. To find this optimal pair of 
models, we use the following procedure: 

1. Optimize true (S) model: FI:(SI— >SC) 

Vary the ionization history of the SI model to find 
the SC model that best matches FI. 

2. Optimize false (F) model: (FI-^FC):SC 

Taking the optimal SC model from step 1 as the 
true model used to generate simulated data, vary 
the ionization history of the FI model to find the 
FC model that matches the best-fit SC. 

Then the significance of rejecting false positives includ- 
ing reionization freedom is —2/S.\nCEE computed for the 
optimal FC:SC pair, i.e. the maximum likelihood from 
step 2. 

These steps for the false positive tests are illustrated 
in Fig. [T2] for the ideal experiment. The process for false 
negatives can be described by simply swapping which 
models have features and which are smooth (F ^ S). 



FIG. 12: False positive example of the two-step process to account 
for reionization uncertainty in polarization significance for the ideal 
experiment. Upper panel: FI:(SI— >SC) — varying the ionization 
history from SI to SC {red curves) to match FI (solid black). Lower 
panel: (FI— >FC):SC — varying the ionization history from FI to FC 
(black curves) to match SC (dashed red). All polarization spectra 
are plotted relative to FI. This joint optimization minimizes the 
FLSI difference [shading in upper panel, with fiducial significance 
given in Table IIIIl l at £ < 30 using the optimal FC:SC models 
(shading in lower panel; Table \V\ . See ^ IIV Bl for an explanation 
of the notation used here. 



Note that even with our conservative choice of Zmax = 
50, Fig. [T^] shows that the main impact of reionization 
on the polarization spectra is limited to ^ < 30. We will 
explore the consequences of this restriction to large scales 
in the following section. 

To implement the steps described above, we vary ion- 
ization histories to minimize —2/S.hiCEE following the 
methods of Refs. H, HO, HO] , using CosmoMC ^ Hi for 
MCMC likelihood analysis with a version of CAMB ^ 
modified to include reionization histories parametrized 
with principal components as in Eq. (j23|) . For example, 
in the FI:(SI^SC) step above, we take the FI model as 
the simulated polarization data and search over ioniza- 
tion histories of the SC model class. 

For each optimization step, we run 4 MCMC chains 
long enough to be well past any initial burn-in phase and 
stop when the region of parameter space near the best fit 
is sufficiently well sampled that all 4 chains agree on the 
maximum likelihood to within ~ 1% in — 2Aln£. Typi- 
cally this requires computing an initial chain to estimate 
the covariance matrix of the reionization PC amplitudes, 
followed by generating chains with ~ 10^ samples each. 
The optimal true or false model is taken to be the overall 
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FIG. 13: Test of false positives due to reionization for an ideal- 
ized experiment limited by sample variance and for Planck. Thick 
curves show the true smooth model {long dashed red) and best- 
fit false feature model {short dashed black) for the false positive 
scenario that would be the most difficult to reject due to free- 
dom in the reionization history. For comparison, the instantaneous 
reionization polarization spectra from Fig. |4] are plotted as thin 
solid curves. Reionization histories are parametrized by 7 principal 
components that cover redshifts 6 < z < 50. 



maximum likelihood model from the final 4 chains. 

All cosmological parameters besides the 7 reioniza- 
tion PC amplitudes are assumed to be fixed by mea- 
surements of the temperature spectrum, except for the 
amplitude of scalar fiuctuations Ag which is varied to 
keep Ag exp(— 2t) fixed, preserving the temperature and 
polarization power at small scales. Fixed parameters are 
set to the values in Table HI We use top-hat priors on the 
PC amplitudes corresponding to < Xe < 1 as described 
in Ref. [i^. Note that the number of PCs used here is 
larger than the 3 to 5 needed for completeness in Ref. [4§| 
due to our choice of a larger maximum redshift. 

Although we are interested in the ability of polarization 
data to test features appearing in the observed temper- 
ature spectrum, we include the contributions from the 
model TT and TE spectra as well as EE in the likeli- 
hood for MCMC. Keeping the temperature data in the 
likelihood ensures that we do not obtain models that fit 
the polarization spectrum well at the expense of chang- 
ing the shape of Cj'^ . For example, ionization histories 
with sharp transitions in Xe at high redshift can gener- 
ate polarization power at £ ~ 40 to match infiationary 
features, but these models also add power to the tem- 
perature on similar scales through an enhanced Doppler 
effect [I^Ig^]. For the best-fit models, the contribution 
of temperature data to the likelihood is approximately 
constant: —2/S.\x\Ctt ~ 7. 



FIG. 14: Same as Fig. 1131 but for tests of false negatives. Here the 
true model has a feature in A^(fc) and the false model is assumed 
to have smooth A|j(fc) (c = 0). 

C. Reionization Confusion 

The optimal true and false model spectra obtained 
from the steps described in the previous section (SC and 
FC) are plotted in Figs. [T51 and [HI for each of our 4 sce- 
narios (ideal/Planck tests of false positives/negatives). 
Wc also plot the corresponding models with instanta- 
neous reionization histories (SI and Fl) to show where 
ionization freedom has the largest effect on the spectra. 
Table|V]lists ~2l\\nCEE for each FC:SC or SC:FC com- 
parison. 

The ability of reionization to cither mimic or obscure 
the signature of inflationary features is greatest in the 
low-power 10 < ^ < 30 regime of Cf^. For tests of 
false positives with an ideal experiment, nearly all of the 
25% reduction of V— 2A In Cee due to reionization comes 
from i < 30. Planck, on the other hand, has relatively 
greater sensitivity to small changes in the polarization 
bump at ~ 40 since observations at such scales suffer 
less from instrumental noise than at i ^ 20. The 17% 
reduction of \/— 2 Ahi Cee for false positive rejection for 
Planck due to reionization is split equally between £ < 30 
and £ > 30. This fact combined with the weakness of 
reionization effects at £ > 30 makes Planck somewhat less 
sensitive to reionization uncertainties than an idealized 
noise- free experiment. 

For false negative tests, changes to polarization spec- 
tra at £ > 30 are generally more important than they are 
for false positives. In the case of the ideal experiment, 
if we ignored multipoles above i = 30 reionization would 
only reduce y/—2A IuCee by 20% relative to the instan- 
taneous reionization significance instead of the 36% re- 
duction that we find when including all scales. The false 
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model being tested in this case (smooth A^(fc)) has less 
power and therefore lower sample variance at 30 < £ < 50 
than the spectrum with a feature, and therefore changes 
in the polarization spectra on these scales have a greater 
effect on the significance than they do for false positive 
tests. Likewise, Planck's significance is more dependent 
on the € ~ 40 bump for testing false negatives than for 
false positives. In fact, nearly all of the 13% degradation 
in \/— 2A IiiCee for false negative rejection comes from 
l>30 for Planck. 

Changes in the reionization history at 6 < 2 < 50 are 
unable to significantly affect the polarization power spec- 
trum at £ > 50. A detection of polarization features on 
these scales would therefore be robust to reionization un- 
certainty. Likewise, measurement of a smooth spectrum 
on these scales would strengthen bounds on the height 
and width of a step in the inflaton potential. 

By considering variations in Xe{z) up to 2 = 50, we 
include a wide variety of ionization histories, many of 
which may not be physically plausible. In practice, how- 
ever, the ionization histories of the spectra in Figs. [T3l 
and [13] have Xe ^ 0.2 at 2 > 20. Nevertheless, had 
we chosen to limit ionization variation to lower redshifts 
the possibility of confusing reionization with inflationary 
features would be lessened, particularly for tests of false 
negatives and for Planck, due to the greater reliance on 
small-scale features in those cases. 

Note that the effects of optimizing the ionization his- 
tory and smoothing the polarization spectra with the 
addition of a large tensor component (§ IIII C|) are simi- 
lar: both are able to make a spectrum with inflationary 
features and a smooth spectrum appear more alike at 
10 ^ ^ ^ 30. Due to this similarity, we expect that con- 
sidering tensors and reionization simultaneously would 
not further degrade the significance of false positive or 
negative tests. 

For variations in the potential parameters discussed 
in § IIIIBI that retain only the dip at £ ^ 20 and not 
the bump a,t £ ^ 40, the impact of reionization will be 
greater. In these cases one cannot expect polarization to 
provide unambiguous confirmation of features without 
external input on the ionization history. 



V. DISCUSSION 

Models with a step in the inflationary potential pro- 
duce oscillations in the angular power spectra of the CMB 
that can improve the fit to WMAP temperature data at 
multipoles £ ^ 20 — 40 at the expense of 3 additional phe- 
nomenological parameters controlling the a step height, 
width, and location on the infiaton potential. Such mod- 
els predict that these oscillations should appear in the E- 
mode polarization spectrum on similar, few-degree scales. 
The first precise measurements of the polarization on 
these scales are anticipated in the next few years from 
Planck, enabling tests of the infiationary-step hypothe- 
sis. 



Moreover, inflationary features at the upper range of 
^ ^ 30 (A: > 2 X 10~3 Mpc'^) that are smoothed out 
due to projection effects in temperature should be more 
visible in polarization. For the lower range of < 30 
{k < 2 X 10^'^ Mpc^^), it becomes important to assess 
the impact of reionization and tensor mode uncertainties. 

We have explored in detail the prospects for polariza- 
tion tests of features, focusing in particular on the risk 
of errors that can be classifled as false positives (falsely 
confirming an infiationary feature) and false negatives 
(falsely rejecting an infiationary feature). Under the sim- 
plest set of assumptions for large-scale polarization in 
which we take the best-fit model for the temperature fea- 
tures, neglect tensor fluctuations, and take the reioniza- 
tion history to be instantaneous, polarization measure- 
ments from Planck should be able to confirm or exclude 
the infiationary features that best match current tem- 
perature data with a significance of ^y—2A IiiCee ~ 3. 
All-sky experiments beyond Planck could potentially in- 
crease this significance to %/~2A IhCee ^ 8, providing 
a definitive test for features from infiation. 

The estimated significance degrades slightly with the 
addition of a large-amplitude, smooth tensor component 
to the i?-mode spectrum, which tends to hide the effect of 
an inflationary step at the largest scales. Assuming that 
the step modifies an to^0^ potential, y/—2A\nCEE is re- 
duced by 4% for Planck and ~ 14% for a cosmic variance 
limited experiment. Allowing non-standard reionization 
histories with arbitrary changes to the ionized fraction 
at 6 < 2 < 50 can lower ^/—2A\nCEE by as much as 
~ 15% for Planck and ~ 30% for cosmic variance limited 
data. Since tensor fiuctuations and reionization have the 
greatest impact on detectability of inflationary features 
at similar scales {£ ~ 20), their effects on the significance 
should not be cumulative. 

The possible contamination due to tensors or reion- 
ization could eventually be mitigated with constraints 
from other types of observations, e.g. stronger limits 
on the tensor-to-scalar ratio from the B-mode polariza- 
tion power spectrum. The -B-mode contribution from 
an vn?( i>^ p otential is potentially within the reach of 
Planck [6J|. Note, however, that a failure to reject false 
positives or negatives for inflationary features in i5-mode 
polarization would generally bias the inferred ionization 
history and reionization parameters such as the optical 
depth. Such biases would in turn lead to biased con- 
straints on inflationary parameters from tensor i?-mode 
measurements [soj . 

These estimated significances assume that the param- 
eters of the step in the infiaton potential are those that 
best fit the WMAP temperature spectrum. Away from 
this best fit, the polarization significance can either in- 
crease or decrease. Cases where the significance substan- 
tially decreases correspond to parameter combinations 
where at most one of the dip {£ ~ 20) and bump {£ ~ 40) 
temperature features can be explained by the step in the 
potential. In the case that only the dip is inflationary, 
Planck will be unable to confirm the feature. 
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We have not computed the impact of foreground re- 
moval uncertainties on our results; in general one might 
expect our forecasts to degrade somewhat upon including 
them. However, recent studies for Planck [64 and a fu- 
ture dedicated polarization satellite mission |65l . f66j indi- 
cate that foregrounds will not be a substantial problem in 
the relevant multipole range. Finally, we do not address 
the possibility that the features in the WMAP data arise 
from a systematic effect (cf. Appendix El . Nonetheless, 
if all of the ^ = 20 — 40 features in the temperature power 
spectrum are inflationary, polarization should ultimately 
provide a statistically significant confirmation. 

Acknowledgments: We thank Jan Hamann and 
Richard Easther for valuable conversations. MJM, CD 
and WH were supported by the KICP through the grant 
NSF PHY-01 14422 and the David and Lucile Packard 
Foundation. MJM was additionally supported through 
the NSF GRFP. WH was additionally supported by the 
DOE through contract DE-FG02-90ER-40560. HVP was 
supported in part by Marie Curie grant MIRG-CT-2007- 
203314 from the European Commission, and by a STFC 
Advanced Fellowship. HVP thanks the Galileo Galilei 
Institute for Theoretical Physics for the hospitality and 
the INFN for partial support during the completion of 
this work. 



APPENDIX A: RELATION TO PRIOR WORK 

Our best-fit parameters for the WMAP temperature 
data (Table |T| differ from those found by previous stud- 
ies of the same inflationary model in Refs. (33l. [34j. We 
explain here the reasons for these discrepancies, which 
include the addition of data and changes in the likeli- 
hood code in going from 3-year to 5-year WMAP data, 
as well as differences in the computation of the evolution 
of modes during inflation. 

Due to small changes in the observed TT spectrum 
between WMAP 3-year and 5-year data, we find that 
the best-fit width of the feature increased from d = 0.022 
to a value of d = 0.027. The lower value of d agrees with 
the best-fit value found by Ref. [131 , which was based on 
the 3-year data. Note that a wider feature in k implies a 
narrower feature in £ for the CMB power spectra. Fig.fTSl 
shows the best fit models for both data sets along with 
the appropriately binned data. 

Updates to the WMAP likelihood code could also cause 
small changes in the best-fit potential parameters. In 
fact, one might be concerned that the feature in the 
WMAP temperature data is only a systematic effect with 
some artificial origin in the likelihood calculation. In par- 
ticular, given the location of a feature, its significance 
could emerge in some fashion from the transition between 



the low-£ pixel-based TT likelihood code and the high-^ 
harmonic space likelihood code, which happens at ^ = 32 
in the 5-year likelihood code However, in the original 
version of the 3-year likelihood code, v2pl, the transition 
occurred at £ = 16, and in the final version, v2p2p2, it 
was changed to £ = 32 [b^]- We searched for the best- fit 
feature model using WMAP3 data with these two ver- 
sions of the likelihood code and found almost exactly the 
same values for the potential parameters in both cases, 
indicating that this particular issue is not the source of 
a systematic effect. 

Our best fit value for b is considerably different from 
Refs. [Hjlll] even though we use the same matching con- 
dition as they do between e-folds and physical wavenum- 
bers. This is due to a choice of initial conditions for the 
background evolution of the inflaton by these authors 
that did not quite satisfy the Friedmann equation, with 
the result that the subsequent evolution also failed to sat- 
isfy it [6^. This essentially translates into a horizontal 
shift in (f>, changing the preferred location of the step b. 

Ref. [131 discusses relaxing the model dependence of 
the predicted power spectrum from the chaotic inflation 
"toy model" adopted here by using a free spectral index 
that is fit to the data rather than set by the choice of iV^,. 
Since the value of w 0.96 determined by our match- 
ing condition for the chaotic inflation potential as de- 
scribed in §|n]is nearly identical to the spectral tilt in the 
WMAP5 best-fit concordance model (i.e. with smooth 
A^(fc)), we do not carry out this extra step here. How- 
ever, the form of the underlying potential will be tested 
by the Planck satellite irrespective of the existence of 
features; as we note in § IIII CI the tensor amplitude pre- 
dicted by the m^0^ potential (which is not affecte d by 
the presence of the feature) is within Planck's reach [64| . 

3000 I I 




FIG. 15: Observed TT spectrum (binned in £, with WMAP3 and 
WMAP5 points offset slightly in £ for clarity) and best-fit feature 
models for WMAP3 (red) and WMAP5 (black). Parameters other 
than c and d are set to the values in Table |I] for WMAP5, and 
to the following values for WMAP3: {m = 6.852 X 10-'',fc = 
14.67, Qth'^ = 0.02222, Qch^ = 0.09927, h = 0.753, r = 0.0817}. 
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